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Abstract 



The goal of the DIRAC experiment at CERN (PS212) is to measure the 7r + 7r™ atom lifetime 
with 10% precision. Such a measurement would yield a precision of 5% on the value of the 
5-wave -kit scattering lengths combination \ao — a 2 \. Based on part of the collected data 



we present a first result on the lifetime, r = 2.91 
systematic errors. This lifetime corresponds to \oq — a 2 
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1 Introduction 



The aim of the DIRAC experiment at CERN [1] is to measure the lifetime of pi- 
onium, an atom consisting of a n + and a n~ meson (A 27V ). The lifetime is dom- 
inated by the charge-exchange scattering process — > 7r°7r°) 1 and is thus 
related to the relevant scattering lengths [4]. The partial decay width of the atomic 
ground state (principal quantum number n — 1, orbital quantum number I — 0) is 
[2,5,6,7,8,9] 

T ls = — = la 3 p\a -a 2 \ 2 (l + 5) (1) 



* PH Division, CERN, CH 1211 Geneva 23, Switzerland; 
e-mail: Leonid.Afanasev@cern.ch 

1 Annihilation into two photons amounts to w 0.3% [2,3] and is neglected here. 
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with Tis the lifetime of the atomic ground state, a the fine-structure constant, p the 
7T° momentum in the atomic rest frame, and a and a 2 the S-wave nix scattering 
lengths for isospin and 2, respectively. The term 5 accounts for QED and QCD 
corrections [6,7,8,9]. It is a known quantity (5 = (5.8 ± 1.2) x 1CT 2 ) ensuring 
a 1% accuracy for Eq. (1) [8]. A measurement of the lifetime therefore allows to 
obtain in a model-independent way the value of \a — a 2 1 . The -kit scattering lengths 
cio, a 2 have been calculated within the framework of Standard Chiral Perturbation 
Theory [10] with a precision better than 2.5% [11] (a = 0.220 ± 0.005, a 2 = 
-0.0444 ± 0.0010, a - a 2 = 0.265 ± 0.004 in units of inverse pion mass) and lead 
to the prediction r is = (2.9 ±0.1) x 10~ 15 s. The Generalized Chiral Perturbation 
Theory though allows for larger a-values [12]. Model independent measurements 
of a have been done using K e4 decays [13,14]. 

Oppositely charged pions emerging from a high energy proton-nucleus colli- 
sion may be either produced directly or stem from strong decays ("short-lived" 
sources) and electromagnetic or weak decays ("long-lived" sources) of interme- 
diate hadrons. Pion pairs from "short-lived" sources undergo Coulomb final state 
interaction and may form atoms. The region of production being small as compared 
to the Bohr radius of the atom and neglecting strong final state interaction, the cross 
section a\ for production of atoms with principal quantum number n is related 
to the inclusive production cross section for pion pairs from "short lived" sources 
without Coulomb correlation (a®) [15] : 



dp A 



= (2tt) 



E 



A hTfC 



dp + dp^ 



M A 



^(f = 0)| 



(2) 



p+=p- 



with pa, Ea and Ma the momentum, energy and mass of the atom in the lab 
frame, respectively, and p + , p_ the momenta of the charged pions. The square of 
the Coulomb atomic wave function for zero distance f* between them in the cm. 
system is |\1/^(0)| 2 = p 3 B /nn 3 , where p B = is the Bohr momentum of the 

pions and m % the pion mass. The production of atoms occurs only in ^-states [15]. 

Final state interaction also transforms the "unphysical" cross section a° s into a real 
one for Coulomb correlated pairs, a c [16,17]: 



dp + dp- 



*_r,(nr^v, o) 



where ^^(P) is the continuum wave function and 2k* = q with q being the 
relative momentum of the 7r + and ir~ in the cm. system 2 . |\I/^(r*)| 2 describes 

2 For the sake of clarity we use the symbol Q for the experimentally reconstructed and q 
for the physical relative momentum. 
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the Coulomb correlation and at r* = coincides with the Gamov-Sommerfeld 
factor A c (q) with q = \q\ [17]: 



A ( ) - 27rm 7r o;/g 

1 — exp(— 2irm 7T a/q) 



For low g, < g < g , Eqs. (2, 3, 4) relate the number of produced A 27T atoms, Na, 
to the number of Coulomb correlated pion pairs, N C c [18] : 



OO 



1 



N A o't (2iram n f ^ n 3 



! j3 g 



fe*h(5o)- (5) 



Eq. (5) defines the theoretical A;-factor. Throughout the paper we will use 

g = 2MeV/c and k th (q ) = 0.615. (6) 



In order to account for the finite size of the pion production region and of the two- 
pion final state strong interaction, the squares of the Coulomb wave functions in 
Eqs. (2) and (3) must be substituted by the square of the complete wave functions, 
averaged over the distance r* and the additional contributions from 7r°7r° — > A 27r as 
well as 7T 7r° — > n + ix~ [17]. It should be noticed that these corrections essentially 
cancel in the A;-factor (Eq. (5)) and lead to a correction of only a fraction of a 
percent. Thus finite size corrections can safely be neglected for k th . 

Once produced, the A 2n atoms propagate with relativistic velocity (average Lorentz 
factor 7 17 in our case) and, before they decay, interact with target atoms, 
whereby they become excited/de-excited or break up. The it + it~ pairs from break- 
up (atomic pairs) exhibit specific kinematical features which allow to identify them 
experimentally [15], namely very low relative momentum q and q^ (the component 
of q parallel to the total momentum p + as shown in Fig. 1. After break-up, the 
atomic pair traverses the target and to some extent loses these features by multiple 
scattering, essentially in the transverse direction, while q L is almost not affected. 
This is one reason for considering distributions in Q L as well as in Q when analyz- 
ing the data. 

Excitation/de-excitation and break-up of the atom are competing with its decay. 
Solving the transport equations with the cross sections for excitation and break-up, 
[20,21,22,23,24,25,26,27,28,29,30,31] leads to a target- specific relation between 
break-up probability and lifetime which is estimated to be accurate at the 1% level 
[22,32,33]. Measuring the break-up probability thus allows to determine the life- 
time of pionium [15]. 
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1 2 3 4 5 

q, \q L \ [MeV/c] 

Fig. 1. Relative momentum distributions (q, q^) for atomic 7r + 7r~ pairs at the point of 
break-up and at the exit of the target. Note that q^ is almost not affected by multiple scat- 
tering in the target. 

The first observation of the A 27r atom [34] has allowed to set a lower limit on its 
lifetime [18,19] of r > 1.8 x 1CT 15 s (90% CL). In this paper we present a deter- 
mination of the lifetime of the A 27r atom, based on a large sample of data taken in 
2001 with Ni targets. 



2 The DIRAC experiment 

The DIRAC experiment uses a magnetic double-arm spectrometer at the CERN 
24 GeV/c extracted proton beam T8. Details on the set-up may be found in [35]. 
Since its start-up, DIRAC has accumulated about 15'000 atomic pairs. The data 
used for this work were taken with two Ni foils, one of 94/im thickness (76% 
of the 7r + 7r~ data), and one of 98 /im thickness (24% of the data). An extensive 
description of the DIRAC set-up, data selection, tracking, Monte Carlo procedures, 
signal extraction and a first high statistics demonstration of the feasibility of the 
lifetime measurement, based on the Ni data of 2001, have been published in [36]. 

The set-up and the definitions of detector acronyms are shown in Fig. 2. The main 
selection criteria and performance parameters [36] are recalled in the following. 

Pairs of oppositely charged pions are selected by means of Cherenkov, preshower 
and muon counters. Through the measurement of the time difference between the 
vertical hodoscope signals of the two arms, time correlated (prompt) events (cr^t = 
185 ps) can be distinguished from accidental events (see [36]). The resolution of the 
three components of the relative momentum Q of two tracks, transverse and parallel 
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MSGC SFD IH 



Magnet VH hh 



Shielding Ch / 

PSh j 
Mu 

Fig. 2. Schematic top view of the DIRAC spectrometer. Upstream of the magnet: target, mi- 
crostrip gas chambers (MSGC), scintillating fiber detectors (SFD), ionization hodoscopes 
(IH) and iron shielding. Downstream of the magnet: drift chambers (DC), vertical and hor- 
izontal scintillation hodoscopes (VH, HH), gas Cherenkov counters (Ch), preshower detec- 
tors (PSh) and, behind the iron absorber, muon detectors (Mu). 

to the cm. flight direction, Q x , Q y and Q L , is about 0.5 MeV/c for Q < 4MeV/c. 
Due to charge combinatorials and inefficiencies of the SFD, the distributions for the 
transverse components have substantial tails, which the longitudinal component 
does not exhibit [37]. This is yet another reason for analyzing both Q and Ql 
distributions. 

Data were analyzed with the help of the DIRAC analysis software package ARI- 
ANE [39]. 

The tracking procedures require the two tracks either to have a common vertex 
in the target plane ("V-tracking") or to originate from the intersect of the beam 
with the target ("T-tracking"). In the following we limit ourselves to quoting results 
obtained with T-tracking. Results obtained with V-tracking do not show significant 
differences, as will be shown later. 

The following cuts and conditions are applied (see [36]): 

• at least one track candidate per arm with a confidence level better than 1 % and a 
distance to the beam spot in the target smaller than 1 .5 cm in x and y; 

• "prompt" events are defined by the time difference of the vertical hodoscopes in 
the two arms of the spectrometer of | At\ < 0.5 ns; 

• "accidental" events are defined by time intervals —15 ns < At < —5 ns and 
7 ns < At < 17 ns, determined by the read-out features of the SFD detector 
(time dependent merging of adjacent hits) and exclusion of correlated n~p pairs. 
[36]; 

• protons in "prompt" events are rejected by time-of-flight in the vertical ho- 
doscopes for momenta of the positive particle below 4 GeV/c. Positive particles 
with higher momenta are rejected; 

• e ± and /x ± are rejected by appropriate cuts on the Cherenkov, the preshower and 
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the muon counter information; 

• cuts in the transverse and longitudinal components of Q are Qt < 4 MeV/c and 
\Ql\ < 15MeV/c. The Qt cut preserves 98% of the atomic signal. The Ql cut 
preserves data outside the signal region for defining the background; 

• only events with at most two preselected hits per SFD plane are accepted. This 
provides the cleanest possible event pattern. 



3 Analysis 

The spectrometer including the target is fully simulated by GEANT-DIRAC [38], 
a GEANT3-based simulation code. The detectors, including read-out, inefficiency, 
noise and digitalization are simulated and implemented in the DIRAC analysis code 
ARIANE [39]. The triggers are fully simulated as well. 

The simulated data sets for different event types can therefore be reconstructed with 
exactly the same procedures and cuts as used for experimental data. 

The different event types are generated according to the underlying physics. 

Atomic pairs: Atoms are generated according to Eq. (2) using measured total mo- 
mentum distributions for short-lived pairs. The atomic 7r + 7r~ pairs are generated 
according to the probabilities and kinematics described by the evolution of the 
atom while propagating through the target and by the break-up process (see [40]). 
These it + it~ pairs, starting from their spatial production point, are then propagated 
through the remaining part of the target and the full spectrometer using GEANT- 
DIRAC. Reconstruction of the track pairs using the fully simulated detectors and 
triggers leads to the atomic pair distribution dn^ c /dQ. 

Coulomb correlated n + n~ pairs (CC-background): The events are generated ac- 
cording to Eqs. (3,4) using measured total momentum distributions for short-lived 
pairs. The generated (/-distributions are assumed to follow phase space modified by 
the Coulomb correlation function (Eq. (4)), dN^/dq oc q 2 x Ac(q). Processing 
them with GEANT-DIRAC and then analyzing them using the full detector and 
trigger simulation leads to the Coulomb correlated distribution dN^Q /dQ. 

Non-correlated n + n~ pairs (NC-background): n + n~ pairs, where at least one 
pion originates from the decay of a "long-lived" source (e.g. electromagnetically 
or weakly decaying mesons or baryons) do not undergo any final state interactions. 
Thus they are generated according to dNff£/dq oc q 2 , using slightly softer momen- 
tum distributions than for short-lived sources (difference obtained from FRITIOF- 
6). The Monte Carlo distribution dNff§ /dQ is obtained as above. 

Accidental -k + it~ pairs (acc-background): it + -k~ pairs, where the two pions orig- 
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inate from two different proton-nucleus interactions, are generated according to 
dN^c/dq oc q 2 , using measured momentum distributions. The Monte Carlo distri- 
bution dN^ /dQ is obtained as above. 

All the Monte Carlo distributions are normalized, ff max (dN^ c / 'dQ)dQ = 
N^ , i = CC, NC, acc, with statistics about 5 to 10 times higher than the 
experimental data; similarly for atomic pairs (n A c ). 



The measured prompt distributions are approximated by appropriate shape func- 
tions. The functions for atomic pairs, F A (Q), and for the backgrounds, F B (Q), 
(analogously for Q L ) are defined as: 



Fa{Q) 



dnf c 



n A ^ dQ srj^ 



rp ((~)\ _ _ ' UV ^ ' GG _ ' iVG ^ iVG 

*BW- N MC dQ + N MC dQ + N MC dQ 



with n A c , Nqq, N^c the reconstructed number of atomic pairs, Coulomb- and 
non-correlated background, respectively, and oo acc the fraction of accidental back- 
ground out of all prompt events N pr . Analyzing the time distribution measured with 
the vertical hodoscopes (see [36]) we find uj acc =l . 1% (7.7%) for the 94 /mi (98 /im) 
data sets [36,37] and keep it fixed when fitting. The x 2 function for Q (analogously 
for Q L ) to minimize is 

X h (^AQ) v + (a A )l + (a B )l 



with AQ the bin width and a A , <y B the statistical errors of the Monte Carlo shape 
functions, which are much smaller than that of the measurement. The fit parameters 
are n A c , Nqq, N^c ( see Eq. (7)). As a constraint the total number of measured 
prompt events is restricted by the condition N pr (1 — uo acc ) = Nqq + N^g + n r A ec . 
The measured distributions as well as the background are shown in Fig. 3 (top). 

The data taken with 94 and 98 /im thick targets were analyzed separately. The total 
number of events in the prompt window is N pr = 471290. 

First, we determine the background composition by minimizing Eq. (8) outside of 
the atomic pair signal region, i.e. for Q > 4MeV/c and Q L > 2MeV/c. For this 
purpose we require n r J c = 0. As a constraint, the background parameters Nqq and 
representing the total number of CC- and iVC-events, have to be the same 
for Q and Ql. Then, with the parameters found, the background is subtracted from 
the measured prompt distribution, resulting in the residual spectra. For the signal 
region, defined by the cuts Q = 4MeV/c and Ql = 2MeV/c, we obtain the total 
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Fig. 3. Top: Experimental Q and Ql distributions after subtraction of the prompt accidental 
background, and fitted Monte Carlo backgrounds (dotted lines). The peak at Q = 4MeV/c 
is due to the cut Qt < 4 MeV/c. Bottom: Residuals after background subtraction. The 
dotted lines represent the expected atomic signal shape. The bin-width is 0.25 MeV/c. 

Table 1 

Fit results (94 and 98 fim targets together, background shapes from Monte Carlo (MC)) for 
the parameters Nq§ (total number of CC-events), N^q (total number of NC events) and 
n^ c (atomic pairs) and deduced results for the number of atomic pairs from the residu- 
als (jf£ sldual } an d the number of CC-background events in the signal region {N^q). MC-a: 
background fit excluding the signal region. MC-b: fit of the entire momentum range includ- 
ing Monte Carlo shape for atomic pairs ("shape fit"). The cuts were at Q cu t = 4MeV/c 
and Ql : cut = 2MeV/c. Q and Qi-distributions were fitted together. The normalized \ 2 
were 0.9 for MC-a and MC-b. 







fjrec 


Arrec 


^.residual 

n A 




N* 


MC-a 


Q 


374022 ± 3969 


56538 


6518 ± 373 




106500 ± 1130 




Ql 


same 


same 


6509 ± 330 




82289 ± 873 


MC-b 


Q 


374282 ± 3561 


56213 




6530 ± 294 


106549 ± 1014 




Ql 


same 


same 




same 


82345 ± 783 



number of atomic pairs, n r X and of Coulomb correlated background events, 
Nq^<. Results of fits for Q and Q L together are shown in Table 1. 



CC-background and NC- or acc-backgrounds are distinguishable due to their dif- 
ferent shapes, most pronounced in the Ql distributions (see Fig. 3, top). Acciden- 
tal and NC-background shapes are almost identical for Q and fully identical for 
Ql (uniform distributions). Thus, the errors in determining the accidental back- 
ground uj acc are absorbed in fitting the NC background. The correlation coefficient 
between CC and NC background is —99%. This strong correlation leads to equal 
errors for N^§ and N^g. The CC-background is determined with a precision better 



9 



than 1%. Note that the difference between all prompt events and the background is 
N pr — Nqq — N^c ~ ^accNpr = 6590, hence very close to the number of residual 
atomic pairs (jf A esidual ) a s expected. This relation is also used as a strict constraint 
for fits outside of the signal region (>), N> - Ngg* - Nfig* - (co acc N pr )> = 
and, hence, the fit requires only one free parameter, . 

Second, the atomic pair signal may be directly obtained by minimizing Eq. (8) over 
the full range and including the Monte Carlo shape distribution F A ("shape fit"). 
The signal strength has to be the same for Q and Q L . The result for the signal 
strength n r A ec as well as the CC-background below the cuts, Nqq, are shown in 
Table 1. The errors are determined by MINOS [41]. 

The consistency between the analysis in Q with the one in Q L establishes the cor- 
rectness of the Q T reconstruction. A 2D fit in the variables (Ql, Qt) confirms the 
results of Table 1 . 



4 Break-up probability 



In order to deduce the break-up probability, P br = n A /N A , the total number of 
atomic pairs tia and the total number of produced A 27T atoms, Na, have to be 
known. None of the two numbers is directly measured. The procedure of obtaining 
the two quantities requires reconstruction efficiencies and is as follows. 

Number of atomic pairs: Using the generator for atomic pairs a large number of 
events, n 9 A en , is generated in a predefined large spatial acceptance window Vt gen , 
propagated through GEANT-DIRAC including the target and reconstructed along 
the standard procedures. The total number of reconstructed Monte Carlo atomic 
pairs below an arbitrary cut in Q, n A c '~ rec (Q < Q cu t) defines the reconstruction 
efficiency for atomic pairs e a A = n A c ~ rec (Q < Q C ut)/n A n . The total number of 
atomic pairs is obtained from the measured pairs by n A = n r A ec (Q < Qcut)/^- 

Number of produced A 27r atoms: Here we use the known relation between pro- 
duced atoms and Coulomb correlated pairs (CC-background) of Eq. (5). Us- 
ing the generator for CC pairs, events, of which N^{q < q ) (see Eq. 6) 
have q below q , are generated into the same acceptance window Vl gen as for atomic 
pairs and processed analogously to the paragraph above to provide the number of 
reconstructed CC-events below the same arbitrary cut in Q as for atomic pairs, 
N cc~ rec (Q < Qcut)- These CC-events are related to the originally generated CC- 
events below g through eg£ = N^~ rec (Q < Q C ut)/N 9 cc(q < <?o)- The number 
of produced atoms thus is Na = hh(qo)N^(Q < Q cu t)l^cc ( see Eq- (6)). 
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Table 2 

k(Qcut) factors as a function of cuts in Q and Q l for the 94 and 98 fim thick Ni targets, and 
the weighted average of the two for a relative abundance of 76% (94fim) and 24% (98/im). 





I7*± JO LLUCI LLL/C 


Q cut = 2MeV/c 


0.5535 ±0.0007 0.5478 ± 0.0007 0.5521 ± 0.0007 


Q cut = 3MeV/c 


0.2565 ± 0.0003 0.2556 ± 0.0003 0.2563 ± 0.0003 


Q cni = 4MeV/c 


0.1384 ±0.0002 0.1383 ±0.0002 0.1384 ± 0.0002 


Q LiCUt = lMeV/c 


0.3054 ± 0.0004 0.3044 ± 0.0003 0.3050 ± 0.0004 


Q L , cut = 2MeV/c 


0.1774 ±0.0002 0.1776 ±0.0002 0.1774 ± 0.0002 



The break-up probability P 6r thus becomes: 

n A n r ^ c (Q < Qcut) , . , „ s. , / \ e< A * /m 

Pl " = n2 = MQ»,)^W < W Wlth ^ = *" (?o) «' <9) 

In Table 2 the /c-factors are listed for different cuts in Q and Q L for the two target 
thicknesses (94fim and 98;um) and the weighted average of the two, corresponding 
to their relative abundances in the Ni data of 2001. The accuracy is of the order of 
one part per thousand and is due to Monte Carlo statistics. 

With the /c-factors of Table 2 and the measurements listed in Table 1, the break-up 
probabilities of Table 3 are obtained. The simultaneous fit of Q and Q L with the 
atomic shape results in a single value. 

The break-up probabilities from Q and Q L agree within a fraction of a percent. The 
values from shape fit and from background fit are in perfect agreement (see Table 
1). We adopt the atomic shape fit value of P br = 0.447 ± 0.023 stat , because the fit 
covers the full Q,Ql range and includes correlations between n^ c and Nqq. 

Analyzing the data with three allowed hit candidates in the SFD search window in- 
stead of two, results in more atomic pairs (see Ref. [36], T-tracking). The break-up 
probabilities obtained are 0.440 ± 0.024 and 0.430 ± 0.021 for Q and Q L , respec- 
tively. They are not in disagreement with the adopted value of 0.447. Despite the 
larger statistics, the accuracy is not improved, due to additional background. This 
background originates from additional real hits in the upstream detectors or from 
electronic noise and cross-talk. This has been simulated and leads essentially to a 
reduced reconstruction efficiency but not to a deterioration of the reconstruction 
quality. The additional sources of systematic uncertainties lead us not to consider 
this strategy of analysis further on. 

V-tracking provides a slightly different data sample, different /c-factors and differ- 
ent signal strengths and CC-background. The break-up probability, however, does 
not change significantly and is p^~ trackin 9 = Q.453 ± 0.025 stat , only 0.3 a off from 
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Table 3 

Break-up probabilities for the combined Ni2001 data, based on the results of Table 1 and 
the /c-factors of Table 2 for the cuts Q cut = 4MeV/c and Qh,cut = 2MeV/c . Errors are 
statistical. 





^residual „rec 

n A n A 


N* 


Pbr 


Q 


6518 ± 373 


106500 ± 1130 


0.442 ± 0.026 


Ql 


6509 ± 330 


82289 ± 873 


0.445 ± 0.023 


Q&Q l 


6530 ± 294 


106549 ± 1004 


0.447 ± 0.023 



the adopted value 0.447. 

The break-up probability has to be corrected for the impurities of the targets. Thus, 
the 94 fim thick target has a purity of only 98.4%, while the 98 fim thick target is 
99.98 % pure. The impurities (C, Mg, Si, S, Fe, Cu) being mostly of smaller atomic 
number than Ni lead (for the weighted average of both targets) to a reduction of the 
break-up probability of 1.1% as compared to pure Ni, assuming a lifetime of 3fs. 
Therefore, the measured break-up probability has to be increased by 0.005 in order 
to correspond to pure Ni. The final result is: 

P br = 0.452 ± 0.023 stat . (10) 



5 Systematic errors 



Systematic errors may occur through the analysis procedures and through physical 
processes which are not perfectly under control. We investigate first procedure- 
induced errors. 

The break-up probability will change, if the ratio N^/N^ depends on the fit 
range. If so, the Monte Carlo distributions do not properly reproduce the mea- 
sured distributions and the amount of CC-background may not be constant. In 
Fig. 4 the dependence is shown for the fits in Q, Q L and both together. The ra- 
tio is reasonably constant within errors, with the smallest errors for a fit range of 
Q = Ql = 15MeV/c. At this point the difference between Q and Q L fits leads to 
a difference in break-up probability of AP^ C = 0.023. 

Consistency of the procedure requires that the break-up probability does not depend 
on Q C ut- In Figure 5 the dependence on the cut is shown for break-up probabilities 
deduced from n r A esidual , There is a systematic effect which, however, levels off for 
large cut momenta. This dependence indicates that the shape of the atomic pair 
signal as obtained from Monte Carlo (and used for the A;-factor determination) is 
not in perfect agreement with the residual shape. This may be due to systematics 



12 



in the atomic pair shape directly and/or in reconstructed CC-background for small 
relative momenta. The more the signal is contained in the cut, the more the Pb r 
values stabilize. As a consequence, we chose a cut that contains the full signal (see 
Eq. (10)). This argument is also true for sharper cuts in Qt than the one from the 
event selection. Cut momenta beyond the maximum cut of Figure 5 would only test 
background, as the signal would not change anymore. 

To investigate whether the atomic pair signal shape is the cause of the above cut 
dependence, we studied two extreme models for atom break-up: break-up only from 
the l.S-state and break-up only from highly excited states. The two extremes result 
in a difference in break-up probability of AP^ ape = 0.008. 




l2rl|MeV/r] 



10 11 12 13 14 15 
Q, \Q L \ [MeV/c] 




3.5 4.0 

e[McV/c] 



Fig. 4. Ratio of CC-background over 
NC -background as a function of fit 
range. 



Fig. 5. as a function of cut momen- 
tum for Q and Ql- 



Sources of systematic errors may also arise from uncertainties in the genuine phys- 
ical process. We have investigated possible uncertainties in multiple scattering as 
simulated by GEANT by changing the scattering angle in the GEANT simulation 
by ±5%. As a result, the break-up probability changes by 0.002 per one percent 
change of multiple scattering angle. In fact we have measured the multiple scatter- 
ing for all scatterers (upstream detectors, vacuum windows, target) and found nar- 
rower angular distributions than expected from the standard GEANT model [42]. 
This, however, may be due also to errors in determining the thickness and material 
composition of the upstream detectors. Based on these studies we conservatively 
attribute a maximum error of +5% and -10% to multiple scattering. 

Another source of uncertainty may be due to the presence of unrecognized K + K~ 
and pp pairs that would fulfill all selection criteria [43]. Such pairs may be as abun- 
dant as 0.5% and 0.15%, respectively, of ix + n~ pairs as estimated for K + K~ with 
FRITIOF-6 3 and for pp from time-of-flight measurements in a narrow momentum 



3 FRITIOF-6 reproduces well production cross sections and momentum distributions for 
24 GeV/c proton interactions. 



13 



Table 4 

Summary of systematic effects on the measurement of the break-up probability P br . Ex- 
treme values have been transformed into a assuming uniform distributions. 



source 


extreme values 


a 


CC -background 


+0.012/ -0.012 


±0.007 


signal shape 


+0.004 / -0.004 


±0.002 


multiple scattering 


+0.01 /-0.02 


+0.006 
-0.013 


K + K~ and pp 


+0 / - 0.04 


+0 

-0.023 


finite size 


+0 / - 0.03 


+0 

-0.017 


Total 




+0.009 
-0.032 



interval with DIRAC data. Their mass renders the Coulomb correlation much more 
peaked at low Q than for pions, which leads to a change in effective it + it~ Coulomb 
background at small Q, thus to a smaller atomic pair signal and therefore to a de- 
crease of break-up probability. The effect leads to a change of AP™' pp = —0.04. 
We do not apply this shift but consider it as a maximum systematic error of P br . Ad- 
mixtures from unrecognized e + e~ pairs from photon conversion do not contribute 
because of their different shapes. 

Finally, the correlation function Eq. (3) used in the analysis is valid for pointlike 
production of pions, correlated only by the Coulomb final state interaction (Eq. 
(4)). However, there are corrections due to finite size and strong interaction [17]. 
These have been studied based on the UrQMD transport code simulations [44] 
and DIRAC data on 7r~7r~ correlations. The parameters of the underlying model 
are statistically fixed with data up to 200MeV/c relative momentum. For Q < 
30MeV/c, the DIRAC data are too scarce to serve as a test of the model. The 
corrections lead to a change of /\pf^ mte " stze _ — 0.02. Due to the uncertainties 
we conservatively consider 1.5 times this change as a maximum error, but do not 
modify P br . 

The systematics are summarized in Table 4. The extreme values represent the 
ranges of the assumed uniform probability density function (u.p.d.f.), which, in 
case of asymmetric errors, were complemented symmetrically for deducing the 
corresponding standard deviations a. Convoluting the five u.p.d.f. results in bell- 
shaped curves very close to a Gaussian, and the ±er (Table 4, total error) corre- 
spond roughly to a 68.5% confidence level and can be added in quadrature to the 
statistical error. 

The final value of the break-up probability is 

P br = 0.452 ± 0.023 stat t° ° f 2 }syst = 0.452 (11) 
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6 Lifetime of Pionium 



The lifetime may be deduced on the basis of the relation between break-up prob- 
ability and lifetime for a pure Ni target (Fig. 6). This relation, estimated to be ac- 
curate at the 1 % level, may itself have uncertainties due to the experimental con- 
ditions. Thus the target thickness is estimated to be correct to better than ± 1 /xm, 
which leads to an error in the lifetime (for P br = 0.45) smaller than ± 0.01 fs, less 
than 1% of the expected lifetime and thus negligible. The result for the lifetime is 

x 10~ 15 s. (12) 



The errors are not symmetric because the Pb r — r relation is not linear, and because 
finite size corrections and heavy particle admixtures lead to possible smaller values 
of Pbr- The accuracy achieved for the lifetime is about +17%, almost entirely due 
to statistics and —21%, due to statistics and systematics in roughly equal parts. 
With full statistics (2.3 times more than analysed here) the statistical errors may be 
reduced accordingly. The two main systematic errors (particle admixtures and finite 
size correction) will be studied in more detail in the future program of DIRAC. 

Using Eq. (1), the above lifetime corresponds to \a — o^l = 0.264 toicio ■ 



p br 

0.48 
0.46 
0.44 
0.42 
0.40 
0.38 
0.36 



1.5 2.0 2.5 3.0 3.5 4.0 

X JS [fs] 

Fig. 6. Break-up probability Pb r as a function of the lifetime of the atomic ground state rig 
for the combined 94 and 98 fim thick Ni targets. The experimentally determined Pbr with 
statistical and total errors translates into a value of the lifetime with corresponding errors. 



ns 



9 qi +0.451 +0.191 

-0.38Js*a* -0A9Jsyst 



X 10 



-15 



2.91 



+0.49 
-0.62 
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